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c3 I We have detected the intrinsic size of Sagittarius A*, the Galac- 

tic Center radio source associated with a supermassive black hole, 
showing that the short-wavelength radio emission arises from very 
near the event horizon of the black hole. Radio observations with 
the Very Long Baseline Array show that the source has a size of 
24 ± 2 Schwarzschild radii at 7 mm wavelength. In one of eight 7- 
mm epochs we also detect an increase in the intrinsic size of 60t^lj%. 
These observations place a lower limit to the mass density of Sgr 
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A* of 1.4 X 10^ solar masses per cubic astronomical unit. 

Sagittarius A* (Sgr A*) is the compact, nonthermal radio source in the Galactic Center 
associated with a compact mass of 4x IO^Mq ()7l[^I^. It is the best established and closest 
supermassive black hole candidate and serves as the prime test case for the black hole 
paradigm. Emission at radio, near-infrared, and X-ray wavelengths traces processes in 
the environment of the event horizon (00001313' • 

High resolution radio imaging of Sgr A* can ultimately distinguish between the many 
different models for the emission, accretion and outflow physics of the source as well as 
provide an important test of strong-field gravity l\10^ . Sgr A* has been a target of such 
observations for the past 30 years ()j Its intrinsic size and structure have remained 
obscured, however, because radio waves from Sgr A* are scattered by turbulent interstellar 
plasma along the line of sight () i^ . The scatter-broadened image of Sgr A* is an ellipse 
with the major axis oriented almost exactly East-West and a quadratic size-wavelength 
relation. 

The turbulent plasma is parametrized with a power-law of turbulent energy density as a 
function of length scale with outer and inner scales that correspond to the scale on which 
turbulence is generated and damped, respectively. Scattering theory predicts that the 
scatter-broadened image will be a Gaussian when the inner length scale of the turbulent 
medium is larger than the longest baseline of the observing interferometer Addition- 
ally, the scatter-broadened image size will scale quadratically as a function of wavelength. 
In the case of Sgr A*, the longest interferometer baseline used in our analysis bmax ~ 2000 
km corresponds to a length scale in the scattering medium D scattering / D 

source ^ ^max ^ 

25 

km, where Dgource = 8 kpc is the distance from Sgr A* to the Earth and Dscattering = 100 
pc is the distance from Sgr A* to the scattering screen This scale is much less 

than the predicted and measured values of the inner scale, which fall in the range 10^ to 
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IQ5.5 j^j^ di^^l The amplitude of turbulence in the Galactic Center scattering screen 
is ~ 2 — 3 orders of magnitude greater than what is seen in the next most powerful scat- 
tering region, NGC 6334B ()j6)). however, suggesting that the Galactic Center case may 
be atypical. 

The presence of strong scattering has pushed observations to shorter and shorter 
wavelengths where scattering effects decrease and intrinsic source structure may dom- 
inate, creating a deviation from the measured size-wavelength law. On the basis of 
extensive observations with the National Radio Astronomy Observatory's Very Long 
Baseline Array (VLBA), L98 measure the index of the size- wavelength power-law to be 
a = 1.99 ± 0.03 ( |j?D . L98 also claim a deviation from the scattering law in the minor 
axis at 7mm wavelength (43 GHz), implying an intrinsic size of 72 Schwarzschild radii 

Unfortunately, precise measurements of the size of Sgr A* are seriously hampered by 
calibration uncertainties related to the variable antenna gain and atmospheric opacity at 
the low antenna elevations necessary to observe Sgr A* from the northern hemisphere. 
Closure amplitudes have been used to constrain the size of Sgr A* with VLBI observations 
at 3.4 mm The closure amplitude does not rely on calibration transfer from another 
source as traditional imaging methods do and is independent of all station-dependent 
amplitude errors. This method does not, however, eliminate baseline-dependent errors 
such as variable decorrelation (which also influence conventional calibration and imaging 
techniques). The closure amplitude is conceptually related to the closure phase, a more 
well-known quantity which is also independent of station-based gain errors. The principle 
drawback of closure amplitude analysis for simple source structures is the reduction in the 
number of degrees of freedom relative to a calibrated data set. The number of independent 
data points for a 7-station VLBA experiment is reduced by a factor 14/21. Additionally, 
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the closure amplitude method can not determine the absolute flux density for the source. 
These shortcomings are more than offset by the confidence that the result gives through 
its accurate handling of amplitude calibration errors. 

We describe here the analysis of new and archival VLBA data through closure am- 
plitude and closure phase quantities. We analyze 3 new experiments including data at 
1.3 cm, 6 new experiments including data at 0.69 cm, as well as 10 experiments from the 
VLBA archive including data at 6, 3.6, 2.0, 1.3, 0.77, 0.69 and 0.67 cm wavelength. 

Observations and Initial Data Reduction 

Six new observations were made with the VLBA as part of our Very Large Array flux 
density monitoring program {BO!)- Three observations were made in each of two sepa^ 
rate epochs in July/August 2001 and April/May 2002 (STable 1). In the first epoch, 
observations at 1.3 cm and 0.69 cm were interleaved over 5 hours. In the second epoch, 
observations were obtained only at 0.69 cm in order to maximize the signal to noise ra- 
tio (SNR) of the final result. All observations were dual circular polarization with 256 
Mbits/sec recording rate. 

We also analyzed a number of experiments from the VLBA archive over the wavelength 
range of 6.0 cm to 0.67 cm (STable 1). The experiments BS055 A, B and C were those 
analyzed by L98. The experiment BB113 was previously analyzed (21). 

Initial data analysis was conducted with the NRAO Astronomical Imaging Processing 
System (f^. Standard fringe-fitting techniques were employed to remove atmospheric 
and instrumental delays from the data (SOM text). High SNR fringes were detected 
for most stations on the compact source NRAO 530 (J1733-1302), indicating the overall 
quality of the data. Due to the relatively larger size of Sgr A*, fringes were obtained for 
a subset of 5 to 8 stations (STable 1). 
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Data were then averaged over wavelength and time for each experiment. The quahty of 
the final result is dependent upon the visibility averaging time. The longer the averaging 
time, the higher the SNR of the closure amplitude calculation (J^^. On the other hand, 
as the averaging time approaches the phase decorrelation time, the closure amplitudes 
cease to be accurate. It is not necessary, however, to determine the best averaging time 
precisely, since neither of these effects is a strong function of time fi£^ . The results that 
we give are for an averaging time of 30 seconds, but we find that for averaging times of 
15 to 120 seconds the estimated intrinsic size of Sgr A* does not differ by more than 10% 
(SOM text). No amplitude calibration was applied at any stage. The averaged data were 
then written to text files for analysis by our own analysis programs, external to AIPS. 

Closure Amplitude and Closure Phase Analysis of a 
Single Gaussian 

We form the closure amplitude from the measured visibilities and average the closure 
amplitudes over time. Closure amplitudes were averaged over scans, which were 5 to 15 
minutes in duration. The code uses the scatter in the closure amplitudes before averaging 
to determine the error in the closure amplitude. Only independent closure amplitudes 
were formed ([^. 

We selected visibility data only with station elevations > 10° to reduce sensitivity to 
phase decorrelation, which is more significant at low elevations. We also excluded data at 
[u, v) distances greater than 25 MX at 6.0 cm, 50 MX at 3.6 cm, 150 MX at 2.0 cm and 
1.3 cm, and 250 MX at 0.69 cm. These sizes are comparable to the expected size of Sgr 
A* at each wavelength. Visibility amplitudes beyond the cutoff were indistinguishable by 
inspection from noise. This {u, f )-distance limit reduced sensitivity to the noise bias or 
station-dependent differences in the noise bias. Results were not strongly dependent on 
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the value of this cutoff. 

Model visibilities for each baseline and time datum were computed for an elliptical 
source of a given flux density 5*0, major axis size x, minor axis size y and position angle 
0. In addition, a noise bias was added in quadrature to each model visibility. Our model 
visibility amplitude (squared) on baseline ij is then 

4 = ^2e-^°(K.-)''-+(-:.i/)^-) + Nl (1) 

where Dn = 2(^^)2, A^",-. is the noise bias, and and f are baseline lengths in 
units of wavelength in a coordinate system rotated to match the position angle 0. Model 
closure amplitudes were then formed from these model visibilities. We determine the best- 
fit parameters using a non-linear fitting method that minimizes between the model 
and measured closure amplitudes (SFig. l,STable 2). We find the reduced for fhe 
amplitudes Xa ~ 1 experiments. 

In the case of an image produced by interstellar electron scattering on baselines longer 
than the inner scale of turbulence, (3 is the power-law index of electron density fluctuations 
(jj,^ . The parameter j3 is related to the exponent a of the scattering law (size oc A°) as 
j3 = a + 2, allowing an independent check of the law (f73|.fmf73j) . For the case of 
the Galactic Center scattering we expect /3 = 4, in which case Equation Q is a Gaussian 
function and x and y are the FWHM in the two axes. Allowing (3 to be unconstrained in 
our fits, we find (3 = 4.00 ± 0.03, which is consistent with the expectation of scattering 
theory (SFig. 3). All remaining analysis is conducted with the assumption that (3 = 4. 

The introduction of the noise bias to the model changes our calculation from a pure 
closure amplitude to a noise-biased closure amplitude. We found that our results did not 
require that we consider the noise bias as dependent on station or time (SOM text). Thus, 
we chose Nij{t) = Nq because it is simpler computationally and has a smaller number of 
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independent parameters. 

Errors in the model parameters were determined by calculating for a grid of models 
surrounding the solution and fitting constant surfaces (SFig. 2). Monte Carlo simula- 
tions find confidence intervals that are smaller by a factor of two than determined from the 

analysis, suggesting that the dominant sources of error are baseline-based errors such 
as phase decorrelation, which were not included in the Monte Carlo simulations (SOM 
text). 

Closure phases were formed, averaged and analyzed in a manner similar to the closure 
amplitudes. We tested the closure phases against the hypothesis that they are all zero. 
This hypothesis is the case for a single elliptical Gaussian and other axisymmetric struc- 
tures with sufficiently smooth brightness distributions. An axisymmetric disk is a notable 
exception to this hypothesis since it induces ringing in the transform plane. The reduced 

for this hypothesis ~ 1 for all experiments (STable 2), indicating no preference for 
multiple components, non-axisymmetric structure or disk-like structure. 

Although the solutions for a single Gaussian component are sufficiently accurate, we 
did search the parameter space for two component models. To do this, we performed 
a minimization of x^ with respect to closure amplitude and closure phase jointly. The 
reduced x^ for these models was roughly equal to the values for the single Gaussian 
component despite the addition of several degrees of freedom. We also calculated upper 
limits to the flux densities of secondary components that are in the range 2-10%, typically 
(SFig. 4, STable 2). The absence of any improvement indicates that a single Gaussian 
component is sufficient and the simplest model of the data. This absence is particularly 
significant for the cases where > 1 and suggests, as noted before, that the results are 
dominated by closure errors rather than improperly modeled structure. 
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Scattering Law and Intrinsic Size 

We determined the size of the major and minor axes of Sgr A* for each experiment (Fig. 1 
and 2, STable 2). The major axis is oriented almost exactly East- West. The major axis 
size is measured much more accurately than the minor axis size because of the poorer 
North-South resolution of the array. All major axis measurements at 1.3 and 0.69 cm are 
larger than the scattering size determined by L98 and the new scattering size that 
we determine below, although the difference is statistically significant in only one epoch 
at 0.69 cm. Minor axis measurements are distributed about the scattering result and no 
one differs significantly from the expected result. 

The L98 scattering law is adequate for the minor axis measurements as a function of 
wavelength (Fig. 3). All the measured minor axis sizes agree with the scattering law to 
better than 3a. The data are also consistent with a constant position angle of 78.0l5;o 
deg with Xu — 2-2 for u = 6 degrees of freedom. 

We determine fits to the major and minor axis sizes as a function of wavelength using 
subsets of the data with a minimum wavelength Xmin of 2.0 cm, 1.3 cm, 0.6 cm and 0.3 
cm (STable 3). The last fit includes the 3.4 mm circular Gaussian fits of for the major 
axis only ()-?.9jl . There are two fits for each subset, allowing a to vary and fixing a = 2. 
xl is less than 3 for the minor axis case with A > 0.6 cm, confirming that the solution is 
adequate for a = 2. 

The major axis data, however, are discrepant from the L98 and the new scattering law 
(Fig. 3). All of the 7mm results fall above the L98 scattering law. Two of these points are 
significantly different at greater than 3cr. The L98 scattering law predicts a size of 690 
/iarcsec at 0.69 cm, which is ~ 7a from the measured size (712I3 yuas) and smaller than 
any of the measured sizes (Fig. 2). An attempt to fit a scattering law with amajor = 2 to 
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all data with A > 0.6 cm gives xl = 24 for 6 degrees of freedom, demonstrating that the 
hypothesis can be strongly rejected. In fact, the 1.35 cm major axis size is also discrepant 
with the best-fit ttmajor = 2 scattering law, giving xl = 5.6 for 3 degrees of freedom. 

We consider two alternative models for our resuls: case A, the scattering power-law 
exponent amajor is not exactly 2; or, case B, intrinsic structure in Sgr A* is distorting the 
size-wavelength relation at short wavelengths. 

For case A, we find adequate solutions for all data at wavelengths > 0.3 cm with 
«major = 1-96 ± 0.01. The result is clearly discrepant with scattering theory which re- 
quires (3 = 4 and marginally discrepant with our determination of the scattering theory 
parameter /3 = 4.00 ± 0.03 (SFig. 3), since scattering theory predicts that a = (3 — 2. 

For case B, we determine a new scattering law from observations with A > 2.0 cm 
and Omajor = 2. This solution has a scale parameter cr^^or ^^at is even less than that of 
L98, increasing the discrepancy at short wavelengths. Removing this new scattering law 
in quadrature gives an intrinsic size of 0.7 ±0.1 mas at 1.35 cm, 0.24 ± 0.01 mas at 0.69 
cm and 0.06 ±0.05 mas at 0.35 cm (Table 1). On the basis of the disagreement between f3 
and a, we reject case A and claim that we have determined the size of intrinsic structure 
in Sgr A* at 1.35 and 0.69 cm. 

The two cases predict substantially different sizes at 20 cm. For the major axis case 
A predicts 541 ± 2 mas while case B predicts 595 ± 3. The 20.7 cm {u = 1450 MHz) 
major axis size 624 ± 6 mas measured with the VLA A-array ()i^6|] is discrepant with both 
of these cases, although more strongly with case A. These measurements are particularly 
difficult since the source is only partially resolved in the A-array: the synthesized beam 
is about 2.6 x 0.9 arcsec oriented North-South. Additionally, extended structure in the 
Galactic Center makes estimation of the size strongly dependent on the estimate of the 
zero-baseline flux density. 
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We attempted to verify the 20 cm size with analysis of three VLA A-array observa- 
tions at 21.6 cm obtained originally for polarimetry (0). Results for each of the three 
experiments were similar and dominated by systematic errors that make an estimate of 
the intrinsic size difficult. We were unsuccessful at analyzing these experiments with our 
closure amplitude technique, possibly due to the poor resolution of Sgr A* and inability 
of our code to handle the large number of stations. In any case, the reliability of ampli- 
tude calibration of the VLA at 20 cm reduces the need for closure amplitude analysis. 
We imaged all baselines and measured the total flux density of Sgr A* by fitting a two- 
dimensional Gaussian to the central 3". For all epochs, we find an error in the total fiux 
density of 10 mJy. We determined the size by fitting in the {u, v) plane with the total fiux 
density fixed and with a minimum cutoff in {u, v) distance. For values of the total fiux 
density that range from —la to +la and for a minimum {u,v) distance from 20 to 120 
kX, we find that the major axis size varies systematically from 580 to 693 mas. The minor 
axis is very poorly constrained. We estimate the size from the mean of these results as 
640 ± 40 mas. We consider this to be a more reasonable estimate of the error in the size 
of Sgr A* than previously given. This size is consistent at < la with case B and ^ 1.5a 
with case A, favoring slightly detection of the intrinsic size. 

Although all minor axis data are adequately fit with aminor = 2, we can check the 
consistency of our results by estimating intrinsic sizes for this axis in the same way. 
The minor axis sizes show the same trend as the major axis sizes: smaller than the 
L98 scattering law at long wavelengths and larger than the L98 scattering law at short 
wavelengths (Fig. 3). Using the solution for aminor = 2 and A > 2.0 cm, we estimate 
intrinsic sizes of 1.1 ± 0.3 mas at 1.35 cm and 0.26 ± 0.06 mas at 0.69 cm. These are 
comparable to the sizes determined for the major axis. For the case of unconstrained 
power-law index fit to all data, we find aminor = l-85lo;o6, marginally consistent with no 
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intrinsic source. 



Changes In the Source Size with Time 

At 0.69 cm, the only measurement deviating significantly from the mean result is in 
the major axis for BB130B. The BB130B resuh is 770l?g /ias while the mean result 
is 712^3 /ias giving a difference of SStig yuas. We note that the greatest deviation in 
the 0.69 cm position angle also occurs for BB130B, although the difference is significant 
only at the 2(7 level. Any such deviation would indicate a non-symmetric expansion or 
a non-symmetric intrinsic source size. We can estimate the change in the size of the 
intrinsic source between BB130B and the mean size by subtracting in quadrature the case 
B scattering size from each. As stated above, the mean result implies an intrinsic size of 
0.24 ± 0.01 mas. The intrinsic size implied by the BB130B result is 0.38to;o4 mas. Thus, 
the growth in major axis size is 0.14l;o;o4 nias in the N-S direction. We cannot associate 
this change in structure with a flux density change. This maximum in the size comes 
~ 10 days before detection of an outburst at 0.69 cm with the VLA {\2(]^ . The following 
epoch, BB130C, occurs only two days before this outburst but shows no deviation from 
the mean size, although the size is particularly poorly determined in this case. 

The Interstellar Scattering Screen 

The image of a scattered source is created by turbulent plasma along the line of sight. 
The minimum time scale for the scattered image to change is the refractive time scale, 
the time in which the relative motions of the observer, turbulent plasma and background 
source lead to the background source being viewed through a completely different region 
of the interstellar plasma. The refractive time scale for Sgr A* is ~ 0.5A^ y cm~^ given a 
relative velocity of 100 km s~^ () At our longest wavelength for VLBA observations, 
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6 cm, then the time scale is 20 y. At our shortest wavelength of 7 mm, the time scale for 
refractive changes is 3 months. Our observations are distributed over a much larger time 
frame than three months, implying that the mean result may be affected by refractive 
changes. 

Two subsets of the archival data have much smaller span, however. The BS055 experi- 
ments cover 6.0 to 0.69 cm in 1 week and the BL070 and BB113 experiments cover 6.0 cm 
to 0.67 cm in 3 months. These data sets include all of the 2.0 cm and longer wavelength 
data. If we compare the 0.7 cm size, we see that it is larger in these quasi-simultaneous 
experiments than in the mean of all experiments and also larger than the expectation of 
the new scattering law. We find 0.69 cm major axis sizes of 728li5 /ias and 713lg^ /ias 
for BS055C and BL070B, respectively, both larger than the mean size of 712^3 /uas (STa- 
ble 2). We conclude that if refractive effects are altering the short wavelength results, 
then their effect is to reduce the deviation from the scattering law, not enhance it. 

Discussion 

Our results allow us to probe the mechanisms responsible for accretion, outflow and 
emission in the vicinity of the black hole. We can compare the measured 7mm intrinsic 
major axis size of 24i?s and its dependence on wavelength with expected values (Fig. 4). 
The intrinsic size of the major axis decreases with wavelength and is best-fit with a power- 
law as a function of wavelength with index Ointrinsic = 1.6 ±0.2. We find for the minor axis 
a similar value aintrinsic = 2.1 ±0.5. Assuming that the source is circularly symmetric and 
using the mean flux density of 1.0 Jy at 7mm 1^20^ . we compute a brightness temperature 



be smaller in the minor axis. A brightness temperature in excess of 10^° K is a strong 
indication that synchrotron radiation is the dominant emission mechanism at work. 




result is a lower limit, because the source may 
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The wavelength-dependent size of Sgr A* now unambiguously shows that the source 
is stratified due to optical depth effects. We rule out models in which the emission 
originates from one or two zones with simple mono-energetic electron distributions l\27^ . 
These models predict a size which is constant with wavelength and is larger than our 
measured size. 

The results are well-fit by a multi-zone or inhomogeneous model, in which the size 
is equal to the radius at which the optical depth is equal to unity {,28}. In a jet model, 
declining magnetic field strength, electron density and electron energy density contribute 
to a size that becomes smaller with wavelength. A detailed jet model for Sgr A* predicts 
an intrinsic size of 0.25 mas at 0.69 cm and 0.6 mas at 1.3 cm (Fig. 4) Exact 
values and wavelength dependence are a function of a number of parameters including 
the relative contributions of the extended jet and the compact nozzle component of the jet. 
The jet model also predicts that the source should be elongated with an axial ratio of 4:1. 
The apparent measured symmetry in the deconvolved sizes in each axis, however, does not 
imply that the intrinsic source is symmetric. For example, an elongated intrinsic source 
that is oriented at 45 degrees to the scattering axis will produce equal deconvolved sizes 
in each axis. Modeling of the closure amplitudes with a complete source and scattering 
model is necessary to determine the elongation for the most general case. 

The thermal, high accretion rate models such as Bondi-Hoyle accretion 1^3 (J^ and ad- 
vection dominated accretion fiows 1^3 l\i require Tg ~ 10^ K, which overpredicts the size in 
each axis by a factor of 3. This disagreement confirms the elimination of these models on 
the basis of the polarization properties of Sgr A* (j2j). On the other hand, the radiatively 
inefficient accretion fiow (RIAF) model has a lower accretion rate and higher Tg, com- 
patible with the polarization and with this measurement. The RIAF model also predicts 
an inhomogeneous electron distribution consistent with a size that reduces with decreasing 
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wavelength. Both the RIAF model and the jet model are similar in the electron energy 
distribution and magnetic field distribution required to produce the observed fiux density 
within the observed size. These models differ principally in the relative contribution of 
thermal electrons to emission in the submillimeter region of the spectrum. 

Extrapolating our size-wavelength relation to longer wavelengths, we estimate a size 
at 2 cm of 130 Rs with a characteristic light travel time of 85 minutes. This is comparable 
to the shortest time scale for radio variability detected, 2 hours, during which the 2.0 cm 
radio flux density changed by 20% (0)- The smooth nature of the spectrum from 90 cm 
to 7 mm, suggests that our size- wavelength relation holds over that entire range 

Our relation implies a size < 2Rs at 1.3mm, comparable to the size of the event 
horizon. The decrease of the source size with wavelength cannot continue much farther 
due to the flnite size of the central object itself. In the millimeter and submillimeter, 
however, the spectral index rises indicating that there may be a break in the size- 
wavelength relation. Ultimately, the size of the event horizon can be viewed as setting a 
limit on the wavelength of the peak emission. The strong break in the spectrum between 
the submillimeter and the NIR may correspond to the wavelength at which the source 
size becomes comparable to the event horizon. Even with a weaker dependence of size 
on wavelength, the light travel time scale at millimeter wavelengths is a few minutes, 
comparable to the shortest time scale observed at X-ray and NIR wavelengths. This 
coincidence suggests that the bright flares observed in at higher energies are 
related to the submillimeter part of the spectrum and come from the vicinity of the black 
hole. The proximity of the millimeter emission indicates that emission at this and shorter 
wavelengths will be subject to strong light bending effects, providing a unique probe of 
strong- fleld general relativity fT^ r ^ . 

The size-wavelength relation also implies that the black hole mass must be contained 



14 



within only a few Schwarzschild radii. Radio proper motion measurements require that 
Sgr A* must contain a significant fraction if not all of the compact dark mass found in 
the Galactic Center (f^[^[^ . Using conservatively only our 7 mm size and the lower 
limit of the Sgr A* mass of 4 x IO^Mq, we find that the mass density in Sgr A* has to be 
strictly above p, > 1.4 x IO^MqAU^'^. The dynamical lifetime of a cluster of objects with 
that density would be less than 1000 years, making Sgr A* the most convincing existing 
case for a massive black hole (T^ . 
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Table 1. Intrinsic Size of the Major Axis of Sgr A 



Wavelength 


Measured Size 


Scattering Size 


Intrinsic Size 


(cm) 


(/ias) 


(/ias) 


(Rs) 


1.35 
0.69 


26351^^ 
712+1 


25331^^ 
669+i 


72+1? 
24+^ 


0.35 


isolig 


1732_2 


6ii 
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2 Measured Size 
Mean of All Experiments 
^ New Scattering Extrapolation 



2.6 



2.7 2.8 2.9 
Major Axis (mas) 



Fig. 1. Sizes from closure amplitude analysis of 3 new and 1 archival 1.3 cm VLBA 
experiments (open circles). The mean size (triangle) is significantly larger than the new 
scattering size (black diamond), which is a fit to all data at A > 2 with a = 2. 
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Fig. 2. Sizes from closure amplitude analysis of 6 new and 2 archival 0.69 cm VLBA 
experiments (open circles). The mean size (triangle) is significantly larger than the new 
scattering size (black diamond). 



21 



CO 

< 

o 

CD 



T3 
N 



1.1 
1.05 



eO.95 
o 

z 0.9 



X J 1- 

< 1.5 



T3 
CD 
N 



0.5 



2 4 
Wavelength (cm) 



16 



1 



0.5 



85 



^80 
< 



8 70 



0.5 



2 4 
Wavelength (cm) 



16 



65 



0.5 



2 4 
Wavelength (cm) 



16 



Fig. S.Major axis size, minor axis size and position angle as a function of wavelength 
normalized to the L98 scattering size. All results are determined by the closure amplitude 
technique except the 21.6 cm result, which is determined from conventional fitting. We 
also include a 3.5 mm measurement in the major axis (fj^. In the upper two panels, solid 
lines show the best-fit scattering model for A > 2 cm in the major and minor axis 
plots. The dotted lines show the best-fit law with a unconstrained for A > 0.6 cm. The 
line in the lower panel is our best fit value of 78 degrees for the position angle. 
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0.1 1 

Wavelength (cm) 

Fig. 4. The intrinsic size of Sgr A* as a function of wavelength. We plot the best-fit size 
in the major axis and 2o upper limits to the size of the minor axis. We also plot one set 
of predictions for the jet length and jet width (f^ . 
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Detection of the Intrinsic Size of Sagittarius A* through Closure 

AmpHtude Imaging: 
Supporting OnUne Material 

Geoffrey C. Bower^, Heino Falclce^, Robeson M. Herrnstein'^'^ Jun-Hui Zhao^, W.M. Goss^ 

& Donald C. Backer^ 

We describe here further details of the analysis of this experiment. 

Fringe fitting analysis in AIPS followed standard practices. Single-band delay solutions 
obtained for NRAO 530 were transferred to Sgr A*. Multi-band delays and rates were 
determined by fringe-fitting Sgr A* itself. Fringes for Sgr A* were found for stations Brewster 
(BR), Fort Davis (FD), Hancock (HN), Kitt Peak (KP), Los Alamos (LA), North Liberty 
(NL), Owens Valley (OV), Pie Town (PT) and St. Croix (SC), although not for all stations 
at all times. For instance, fringes with BR are detected at all wavelengths but only in the case 
of substantial fore-shortening at the longer wavelengths. While fringes were found between 
SC and HN, fringes were not found between either of these stations and any other station 
forcing us to exclude SC and HN from further analysis. We list all stations included in the 
analysis in STable 1. 

We compute the closure amplitude for four stations m, n,p and q: 



where is the amplitude of the visibility on the baseline between stations i and j. The 
closure amplitude is closely related to the closure phase for three stations m, n and p: 



■'mnp 



+ 0np + 0pm- (2) 
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Here (f)ij is the visibility phase on the basehne between stations i and j. 

The noise bias principally has the effect of increasing model visibility amplitudes on long 
basehnes where the Gaussian source is heavily resolved. Our method of including the noise 
in the model visibilities avoids the problems of unbiasing the measured visibility amplitudes 
or the measured closure amphtudes. The noise in the visibihty amphtudes biases the closure 
amplitudes by a factor that depends on the unbiased visibility amplitude on that baseline 
(?). The factor is accurate only for high SNR data and, therefore, not applicable for our 
case because our data is in the domain where SNR~ 1 — 10. A consequence of our technique 
is that we are computing a quantity that differs slightly from the true closure amplitude. 

Time- and station-dependence of the noise bias is unimportant. As the SNR approaches 
1, the closure amplitudes have the largest error and therefore the least weight in the x^-fitting. 
In the case that the noise bias is station- and time-independent, then the result is dependent 
solely on the ratio R = So/N, where 5*0 is the total flux density and N is the noise bias for 
a given integration time, typically 15 minutes. This ratio ranged from i? ~ 10 at 0.69 cm to 
it! ~ 100 at 6.0 cm. The flux density and the noise bias are not physical values since we have 
not performed any amplitude calibration. We found from examination of the data and Monte 
Carlo simulations that the ratio is determined with about 5% accuracy. More importantly, 
the major axis, minor axis and position angle error estimates are not strongly dependent on 
R; that is, these parameters are not covariant with R. For example, the best-fit major axis 
value for the experiment BB130D is TOSlj^g /las. At the 3a extrema in R, the best-fit major 
axis values are 697 and 704 /las, only marginally different. This lack of covariance permits 
us to simplify our error methodology by dropping R and considering only x, y and (p. 

For the final results, we calculated 101^ grid points in major axis, minor axis and position 
angle with a resolution of 5 /larcsec in major axis, 8 /larcsec in minor axis and 0.9 degree in 
position angle at 0.69 cm. Grid spacings in major and minor axis increased with the square 
of wavelength for the other data sets. Significance levels were determined assuming Gaussian 
statistics and computing the appropriate increase in for the number of degrees of freedom 

(?)■ 

Our method for the determination of errors is not strictly correct since the error distri- 
butions for the parameters are not normal. This deviation can be seen in the shces through 
the surface in SFigure 2. If the distributions were normal, then the 3-dimensional error 
surface would be an ellipsoid. Errors for each parameter are determined from the projection 
of the 3-dimensional surface into one dimension. We report errors based on the 99.73% 
confidence interval which corresponds to Scr in the normal case (STable 2). When we cite 
errors throughout the rest of the paper, we define the la error as one-third of the 3a error, 
not the la error determined from the surface. We do this because the la error is not as 
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well determined from the surface. 

We performed Monte Carlo simulations to relate Ax^ to a confidence interval for the 
0.69 cm experiment BB130D. We modeled errors in the data as Gaussian errors in the 
visibility amplitude after vector averaging in wavelength and time. We use the best-fit noise 
to determine the noise distribution. We also include station-dependent gain fluctuations 
as well as a time- and station-dependent gain that represents opacity errors for a uniform 
atmosphere. In the limit that these errors do not change in closure amplitude averaging time 
(e.g., 15 minutes), the closure amplitude method should be independent of these changes. We 
used the best Gaussian parameters for the BB130D data (STablc 2) as an input model for the 
visibilities and performed 1000 iterations of the Monte Carlo test. The resulting parameters 
fit to the data were distributed in a Gaussian fashion and give means of 704 ± 7 //as in the 
major axis, 361 ± 22 //as in the minor axis and 81.4 ± 1.3 deg in the position angle {la 
errors). The mean values are all consistent with the input model but the errors are lower 
by a factor of two to five from the errors determined by fitting. These Monte Carlo 
results set a firm lower limit on the errors. The apparent Gaussian nature of these errors 
and their small values relative to the errors determined from the data, however, suggest 
that our Monte Carlo method does not fully account for sources of error in the data. Phase 
decorrelation effects not included in the Monte Carlo simulation are the most likely source 
of the additional error. 

All sohitions discussed in the paper assume = 4, as required by scattering theory. We 
did, however, perform fits with j3 unconstrained. Not surprisingly, the other parameters were 
determined with less accuracy when we added the additional parameter. Nevertheless, their 
values were comparable with those previously determined. These individual measurements 
of P had an error of 0.05 to 0.2. 

We determine the mean value of /9 as a function of wavelength (SFigure 3). The mean 
for all experiments is j3 = 4.00 it 0.03, fully consistent with the assumption of Gaussian 
nature. There is no evidence for a trend in /3 with wavelength. We conclude that the images 
are precisely fit by a Gaussian intensity distribution, constraining the scattered image to 
strictly follow a A^-law. Any deviation from a A^-law must indicate intrinsic structure. 

Although the solutions for a single Gaussian component are sufficiently accurate, we 
did search the parameter space for two component models. This search requires joint min- 
imization of closure phase and closure amplitude. In all of the models considered, the two 
Gaussians are assumed to be identical in major axis, minor axis and position angle, which is 
the expectation of the scattering model for the Galactic Center on these small angular scales 
(?). For the most general case of this kind involving the relative position, where the three 
size parameters and the fiux densities of the components are free parameters, the solution 
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from our minimization routine was highly sensitive to the initial guess. Accordingly, we 
considered models in which the size and position angle of the components are fixed at the 
L98 scattering value and the relative position of the second component is fixed with respect 
to the main component ?. At 0.69 cm, the scattering model gives a size of 690 x 366 /las 
in a position angle of 80 degrees. Only the two fiux density parameters and the noise bias 
parameter were variable for each model. A 51 x 51 grid was constructed in relative position 
with steps equal to half the scattering size in the given dimension. In the case of 0.69 cm 
images, this corresponds to a field of 17.6 x 8.8 mas. This modeling gives the best-fit fiux 
density of a second component as a function of position. 

In SFigure 4, we show four examples of flux density ratio maps. These results can be 
considered an upper limit on any second component. For the case of BB130A, which is 
typical of the 0.69 cm experiments, the fiux density ratio peaks within the scattering size at 
~ 10%. The apparent double structure is spurious since the fiux density ratio for the central 
point is undefined. The interpretation of this plot is that we cannot discriminate between 
the single component Gaussian and two Gaussians separated by less than half the scattering 
size with a fiux density ratio of 10% or less. 

We list in STable 2 the maximum values of the fiux density ratio, F2, in the maps, 
excluding the central few pixels that arc interior to the scattering region of the first source. 
We find that the maximum ratio at a given wavelength strongly scales with wavelength. The 
limits on secondary components are typically weaker at long wavelengths, which is almost 
certainly a function of the the more limited (u, f )-covcragc available for these experiments. 
At 0.69 cm, limits on a secondary component are typically 5% while at 1.3 cm the limits are 
on the order of 10%. 

Due to the potentially significant effects of phase de-correlation, we performed our com- 
plete analysis for averaging times of 15, 30, 60 and 120 seconds, providing a test of the 
dependence of the result on phase decorrelation. We find sizes of 22 ± 1, 23 ± 2, 24 ± 2 
and 24 ± 2 i?s for the intrinsic size of Sgr A* at 0.69 cm for each of these averaging times, 
respectively. Thus, the results are not biased by our averaging time at a level more than 
IRs- 

We also show that no one experiment determines the significance of the results. Drop- 
ping experiment BB130F, which has the smallest errors, we find that the mean 0.69 cm size 
is 7I6IJ2 A*sec. This size is only slightly larger than the mean size of 712l;g^ //sec. The re- 
sulting discrepancy with the L98 scattering law is 6a instead of 7a. Dropping both BB130F 
and BB130B from the average, we find a size 714l]^2 A*sec, which is also strongly inconsistent 
with the L98 scattering law. 



- 5 - 



We summarize in STable 3 the results of fitting the major and minor axis sizes as a 
function of wavelength. We solve for the power-law index, a, and the normalized size at 
a wavelength of 1 cm, a^"'^, for each axis as a function of the minimum wavelength, Xmin, 
included in the fit. Fits are included with a constrained to 2 and unconstrained. We also 
compute the reduced xl ior each fit. 
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Fig. 1. — Closure amplitude for four station groups from the 0.69 cm experiment BB130D. 

The numbers after the station identification are the reduced the L98 scattering model 

(dashed line) and the best fit model (solid line). The improvement in the model fit is readily 
apparent for these baselines. The station group BR-FD-OV-PT appears twice because two 
independent closure amplitudes can be formed from the set of four stations. These have 
different extent in time because some baselines exceed the 2 SOMA cutoff for baseline length. 
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Phi= 81.2 Minor Axis= 0.364 




Minor Axis (mas) Axial Ratio 

Fig. 2. — Three cuts through the surface in major axis, minor axis and position angle 
space for 68%, 95.4%, and 99.73% confidence intervals of source parameters for one 0.69 
cm observation (BB130D) analysed with the closure amphtude method. These contours 
correspond to 1, 2 and 3a for the case of Gaussian errors. Cuts are centered on the best-fit 
values of each parameter, which are listed above the plot. The lower right corner shows grid 
points of 99.73% confidence interval in axial ratio versus position angle space. 
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Wavelength (cm) 



Fig. 3. — The power-law index of density fluctuations, (3, plotted as a function of wavelength. 
Individual results and the mean value over all wavelengths < /? >= 4.00 ±0.03 are consistent 
with strong scattering on baselines shorter than the inner scale of turbulence. 
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Fig. 4. — Maps of the best-fit flux density ratio for a two Gaussian component model as a 
function of relative position. Contour intervals are 1, 2, 4, 8 and 16%. The maps are for the 
0.69 cm experiments BB130D (upper left) and BB130F (upper right), 1.3 cm experiment 
BB130A (lower left) and 2.0 cm experiment BS055C (lower right). 
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Table 1. Observations of Sgr A* 



Code 


Date 
(dd/mm/yy) 


Freq. 
(GHz) 


BR 


FD 


KP 


LA 


NL 


OV 


PT 


Yl 




New Observations 


BB130A 


12/07/01 


43 






















22 


















BB130B 


29/07/01 


43 






















22 


x/ 












x/ 




BB130C 


05/08/01 


43 


v/ 




















22 


















BB130D 


15/04/02 


43 


















BB130F 


03/05/02 


43 


















BB130G 


13/05/02 


43 


%/ 


%/ 


%/ 














ArcliiA'al Datf 




















BS055A 


07/02/97 


4.98 






















8.42 


















BS055B 


12/02/97 
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15 


















BS055C 


14/02/97 


43 


















BL070B 


23/05/99 


43 


















BL070C 


31/05/99 


45 










x/ 


x/ 










39 








v/ 


v/ 


v/ 






BB113 


29/08/99 


8.42 










v/ 


v/ 










4.99 










v/ 


v/ 
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Table 2. Closure Amplitude and Closure Phase Modeling Results 



Code 


V 

(GHz) 


X y 
(lias) (/^as) 


(deg) 


Xa 


xl 


F2 




5.0 GHz 


BS055A 


5.0 


48638tli 153811^^1? 


87.11i^| 


2.0 


1.7 


0.31 


BB113 
ALL 


5.0 
5.0 




7Q+28 
'"^-15 
oo 9+10.2 


2.2 


1.0 


0.02 


8.4 GHz 


BS055A 


8.4 


17697+^°? 7804+^^^^ 


81.8+1-5 


0.9 


3.3 


0.03 


BB113 


8.4 


16840l?^°o 10150™ 


831}^ 


2.1 


1.3 


0.17 


ALL 


8.4 


1762611^ 8514^1^^^ 


81.9lli 








15 GHz 


BS055B 


15.3 


OiyO_228 -'-yt)0_^9g3 


84.0l}}:2^ 


1.4 


1.1 


0.03 


22 GHz 


BB130A 
BB130B 


22.2 
22.2 


981 7+683 -1 700+468 
31T 1732 

26551?^^ 17361^^^ 


70.91^^:1 

74.2+^f^ 


2.7 
2.0 


1.3 
0.9 


0.11 
0.09 


BB130C 


22.2 


28521^}^ 1721^^1 


es.olls^:? 


2.2 


1.2 


0.15 


BS055B 
ALL 


22.2 
22.2 


2607+^^3 1164+f^2^ 
26351^^2' 1613111^ 


77.8li6o:I 
74.61 


0.6 


1.8 


0.01 


39 GHz 


BL070C 


39.1 


8841 6101^1° 




1.7 


1.0 


0.03 


43 GHz 


BB130A 


43.2 


725+ 11 m+lll 


7P- 7+13.3 
' ^- ' -26.3 


1.3 


1.1 


0.03 


BB130B 


43.2 


7701 11 465l^«° 


69.61^^:° 


1.0 


1.1 


0.06 


BB130C 


43.2 


70411^^ 323li3^ 




1.7 


0.9 


0.07 


BB130D 


43.2 


7081 11 364111 


OI 9+12.3 
Oi.A_i2 9 


0.8 


1.3 


0.04 


BB130F 


43.2 


7081 11 SQOtlTs 


81.41 


0.9 


1.2 


0.03 


BB130G 


43.2 


7091 11 3501111 


81.41 li 


0.7 


1.3 


0.03 


BL070B 


43.1 


7131 11 4341^29 


76.6li|| 


0.9 


1.0 


0.02 


BS055C 


43.2 


70fi+ 47 /IOfi+157 
1 33 4^«_428 


75.8li^7-i 


0.5 


1.9 


0.03 


ALL 


43.2 


7121 407li^^ 


79.81 li 








45 GHz 


BL070C5 


45.1 


683l^f 4581^^^ 




1.2 


1.3 


0.03 



Table 2 — Continued 



Code V X y <!> Xa xl F2 

(GHz) (//as) (//as) (deg) 



Note. — a; is the major axis size of the single elliptical Gaussian component; y is the 
minor axis size; is the position angle; Xa is the reduced ioi the closure amplitudes for 
a single component elliptical Gaussian; X4, is the reduced x^ for the closure phases for an 
axisymmetric model; and, F2 is the maximum ratio of the flux densities in a two component 
model. Errors are 3a. 



Table 3. Fits to the Size of Sgr A* as a Function of Wavelength 



-^min Q^major '-'^major Xi/ d.O.f. CKminor '-'^minor Xi/ d.O.f. 

(cm) (mas) (mas) 



2.0 


2 01° °^ 


1 070.05 

^■"^'0.04 


1.6 


1 


2 01°-i^ 
'^■'-'-'-0.22 


640-25 

'-'■"^0.08 


1.1 


1 




2 




1.0 


2 


2 


0.650:0| 


0.5 


2 


1.3 






2.9 


2 


' ^0.09 


0.910:08 


1.1 


2 




2 


1 40°-°^ 


5.6 


3 


2 


0.680:03 


2.0 


3 


0.6 


-■-■^"o.oo 


1 4fi°-°i 
J--^"o.oo 


2.5 


5 


-■-•0^0.06 


n 81 0.03 

U-8io.o5 


0.9 


5 




2 


1 440.01 

J- -^^0.00 


18.3 


6 


2 


710 03 

^- ' -'-0.04 


2.3 


6 


0.3 


-•-•^^0.00 

2 


1 4^50.01 
-•-•^"o.oi 

1 440.01 

-■-•^^0.01 


2.1 
15.7 


6 
7 











